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We discuss charm production at LHC. The production of single cc pairs is calculated in the k t - 
factorization approach. We use several unintegrated gluon distributions from the literature. The 
hadronization is included with the help of fragmentation functions found for the production of 
charmed mesons in e + e~ collisions. Differential distributions for several charmed mesons are 
presented and compared to recent results of the ALICE and LHCb collaborations. Furthermore 
we discuss production of two pairs of cc within a simple formalism of double-parton scattering 
(DPS). Surprisingly large cross sections, comparable to single-parton scattering (SPS) contribu- 
tion, are predicted for LHC energies. Both total inclusive cross section as a function of energy and 
differential distributions are shown. We include recently discussed evolution of double partons in 
the case of two scales. We discuss perspectives how to identify the double scattering contribution. 
We find much larger cross section for large rapidity distance between charm quarks from different 
hard parton scatterings compared to single scattering. 
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1. Transverse momentum spectra of open charm mesons at LHC 

Recently ALICE and LHCb collaborations have measured inclusive transverse momentum 
spectra of open charm mesons in proton-proton collisions at y/s = 7 TeV [[IJ, W. These measure- 
ments are very interesting from the theoretical point of view because of the collision energy never 
achieved before and unique rapidity acceptance of the detectors. Especially, results from forward 
rapidity region 2 < y < 4, obtained by LHCb can improved our understanding of pQCD production 
of heavy quarks. 

The inclusive production of heavy quark/antiquark pairs can be calculated in the framework 
of the k t -factorization [J3|] . In this approach transverse momenta of initial partons are included 
and emission of gluons is encoded in a so-called unintegrated gluon, in general parton, distribu- 
tions (UGDFs). In the leading-order approximation (LO) within the ^-factorization approach the 
differential cross section for the QQ can be written as: 

da ^ fd 2 K U d 2 K 2 , t T-ZTvi /j Y) 



dy\dp\ t dy2dp2 t d§ j~j J % % \6k 2 {x\X2s) 2 



',/ 

8 2 (it 1)t + K 2 ,t - Pl.t - Pl.t) &i{x\ , K?,) &j(x 2 , K2 >t ) 



where &i(x\, K 2 t ) and JPjfa, K% t ) are the unintegrated gluon (parton) distribution functions. 

There are two types of the LO 2 — > 2 subprocesses which contribute to heavy quarks produc- 
tion, gg — > QQ and qq — > QQ. The first mechanism dominates at large energies and the second one 
near the threshold. Only gg — > QQ mechanism is included here. We use off-shell matrix elements 
corresponding to off-shell kinematics so hard amplitude depends on transverse momenta (virtuali- 
ties of initial gluons). At very high energies, especially in the case of charm production at forward 
rapidities, rather small ^-values become relevant. In this kinematical regime calculation of unin- 
tegrated parton distributions is not under full theoretical control and can include in different ways 
various theoretical aspects, like effect of small-* saturation or treatment of nonperturbative region. 
In order to show the uncertainty of our predictions resulting from different approaches in calculat- 
ing unintegrated parton distributions we have used several models from the literature. All of them 
have different theoretical background. It is therefore very interesting to compare such results with 
the recent ALICE and LHCb data and verify applicability of these UGDFs at LHC energies. More 
details of theoretical model applied here can be found in Ref. [ffl]. 

The hadronization of heavy quarks is usually done with the help of fragmentation functions. 
The inclusive distributions of hadrons can be obtained through a convolution of inclusive distribu- 
tions of heavy quarks/antiquarks and Q — > h fragmentation functions: 



do(y h ,pt. h ) _ r l dz D .. <te^g(yQ 1 Pt& 
dyhd 2 Pt.h Jo z 2 dy Q d 2 p t .Q 



(1.2) 



yQ=y>, 



where p t Q = where z is the fraction of longitudinal momentum of heavy quark earned by 
meson. We have made approximation assuming that yQ is unchanged in the fragmentation process. 

In Fig. |] we present our predictions for differential distributions in transverse momentum of 
open charm mesons together with the ALICE (left panel) and LHCb (right panel) experimental 
data. We plot here results obtained with different models of UGDFs. The calculations of charm 
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quarks are performed for m c = 1.5 GeV with the values of the renormalization and factorization 
scales taken to be /i| = [ip = mj . We get very good description of the experimental data, in both 
ALICE and LHCb cases only with KMR model of unintegrated gluon distributions. The other 
of the applied parametrizations of UGDFs do not work and clearly underestimate experimental 
data points. One can also observe huge difference between results of standard collinear LO partem 
model predictions (dotted line) and those obtained in the LO & r -factorization approach (solid line). 

In FigFJ we discuss uncertainties due to the modification of the charm quark mass m c £ (1.2, 
1.8 GeV) (left panel), as well as uncertainties related to different models of heavy quark fragmen- 
tation (right panel). As one can observe, there is only a small sensitivity of the results to the value 
of charm quark mass at small transverse momenta. Some small uncertainties due to the choice of 
fragmentation function appear only at larger values of meson p t 's. 
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Figure 1: Transverse momentum distributions of D mesons together with the ALICE (left) and LHCb (right) 
data. Each curve corresponds to different model of UGDFs. The calculations were done with the help of 
Peterson fragmentation function. 
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Figure 2: Uncertainties of our predictions related to the charm quark mass (left) and to the choice of frag 
mentation function (right). 



2. Double charm production via Double Parton Scattering 

The mechanism of double-parton scattering (DPS) production of two pairs of heavy quark and 
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heavy antiquark is shown in Fig. || together with corresponding mechanism of single-scattering 
production. The double-parton scattering has been recognized and discussed already in seventies 
and eighties. The activity stopped when it was realized that their contribution at center-of-mass 
energies available then was negligible. Nowadays, the theory of the double-parton scattering is 
quickly developing (see e.g. ||, |7|]) which is partly driven by new results from the LHC. 

In the present analysis we wish to concentrate on the production of (cc) (cc) four-parton final 
state which has not been carefully discussed so far, but, as will be shown here, is particularly 
interesting especially in the context of experiments being earned out at LHC and/or high-energy 
atmospheric and cosmogenic neutrinos (antineutrinos). 

The double-parton scattering formalism in the simplest form assumes two single-parton scat- 
terings. Then in a simple probabilistic picture the cross section for double-parton scattering can be 
written as: 

<T DPS (pp -> ccccX ) = —— o SPS (pp-> ccX { )-e SPS (pp-> ccX 2 ) . (2.1) 

2 °eff 

This formula assumes that the two subprocesses are not correlated and do not interfere. At low 
energies one has to include parton momentum conservation i.e. extra limitations: x\ + x^ < 1 and 
X2+X4 < 1, where x\ and X3 are longitudinal momentum fractions of gluons emitted from one 
proton and X2 and X4 their counterparts for gluons emitted from the second proton. The "second" 
emission must take into account that some momentum was used up in the "first" parton collision. 
This effect is important at large quark or antiquark rapidities. Experimental data M] provide an 



estimate of a e /f in the denominator of formula (|2JJ). In our analysis we take a e /f = 15 mb. 




Figure 3: SPS (left) and DPS (right) mechanisms of (cc)(cc) production. 

A more general formula for the cross section can be written formally in terms of double-parton 
distributions (dPDF), e.g. F gg , F qq , etc. In the case of heavy quark production at high energies: 



d(y DPS = -^— F ' m (xuX3,tf>V$)Fgg(x2,X4,fJLf,fjg] 

2 °eff 



da gg ^ C c{x l ,x 2 ,n l )da gg ^ C c(x3,x 4 ,ij£) dxidx 2 dx 3 dx 4 . (2.2) 

It is physically motivated to write the double-parton distributions rather in the impact parame- 
ter space F gg (xi,X2,b) = g(x\)g(x2)F (b), where g are usual conventional parton distributions and 
F(b) is an overlap of the matter distribution in the transverse plane where b is a distance between 
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both gluons [Q. The effective cross section in (2.1) is then \ja e ff = Jd 2 bF 2 (b) and in this ap- 
proximation is energy independent. 

The double-parton distributions in Eq.(2.2) are generally unknown. Usually one assumes a 
factorized form and expresses them via standard distributions for SPS. Even if factorization is valid 
at some scale, QCD evolution may lead to a factorization breaking. For some time the evolution 
was known only when the scale of both scatterings is the same [0] i.e. for heavy object, like double 
gauge boson production. Recently the evolution of dPDF was discussed also in the case of different 
scales [0]. 

In the left panel of Fig. ^ we compare cross sections for the single cc pair production as well 
as for single-parton and double-parton scattering cccc production as a function of proton-proton 
center-of-mass energy. At low energies the conventional single cc pair production dominates. The 
cross section for SPS production of cccc system [|ll|] is more than two orders of magnitude smaller 
than that for single cc production. For reference we show the proton-proton total cross section as 
a function of energy as parametrizes in Ref. [12]. At low energy the cc or cccc cross sections are 
much smaller than the total cross section. At higher energies the contributions approach the total 
cross section. This shows that inclusion of unitarity effect and/or saturation of parton distributions 
may be necessary. At LHC energies the cross section for both terms becomes comparable. This is 
a new situation when the DPS gives a huge contribution to inclusive charm production. 




Vs (GeV) p. (GeV) 



Figure 4: Total LO cross section for single cc pair and SPS and DPS cccc production as a function of center- 
of-mass energy (left panel) and transverse momentum (right panel) of c or c quarks at yfs = 7 TeV. Cross 
section for DPS should be multiplied in addition by a factor 2 in the case when all c (c) are counted. We 
show in addition a parametrization of the total cross section in the left panel. 



In the right panel of Fig. |] we present transverse momentum distributions of single c (c). 
Within approximations made in this paper the distributions are identical in shape to single-pair 
production distributions. This means that the double-scattering contribution produces naturally 
an extra center-of-mass energy dependent K factor to be contrasted with approximately energy- 
independent ^-factor due to next-to-leading order QCD corrections. Other interesting conclusions 
can be obtained by studying correlation observables. 

In Fig. |5] we show distribution in the difference of c and c rapidities ydiff =y c —yc (left panel) 
as well as in the cc invariant mass M cS (right panel). We show both terms: when cc are emitted in the 
same parton scattering (ciq or C3C4) and when they are emitted from different parton scatterings 
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(C1C4 or C2C3). In the latter case we observe a long tail for large rapidity difference as well as at 
large invariant masses of cc. 

In particular, cc (or cc) should be predominantly produced from two different parton scatter- 
ings which opens a possibility to study the double scattering processes. A good signature of the 
cccc final state is a production of two mesons, both containing c quark or two mesons both contain- 
ing c antiquark (D D or/and D°D°) in one physical event. More detailed discussion of the DPS 



charm production can be found in our original paper Ref. [13]. 
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Figure 5: Distribution in rapidity difference (left panel) and in invariant mass of the cc pair (right panel) at 
V?=7TeV. 

In the present approach we have calculated cross section in a simple collinear leading-order 
approach. A better approximation would be to include multiple gluon emissions. This can be done 
e.g. in soft gluon resummation or in £ r -factorization approach. This will be discussed in detail 
elsewhere. 



References 

[1] B. Abelev et al. (The ALICE Collaboration), J. High Energy Phys. 01 (2012) 128. 

[2] The LHCb Collaboration, LHCb-CONF-20 10-013. 

[3] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. 366 (1991) 135. 

[4] M. Luszczak, R. Maciula and A. Szczurek, Phys. Rev. D79 (2009) 034009. 

[5] A.M. Snigirev, Phys. Rev. D68 (2003) 114012. 

[6] V.L. Korotkikh and A.M. Snigirev, Phys. Lett. B594 (2004) 171. 

[7] J.R. Gaunt and W.J. Stirling, J. High Energy Phys. 03 (2010) 005. 

[8] F. Abe et al. (The CDF Collaboration), Phys. Rev. D56 (1997) 381 1. 

[9] G. Calucci and D. Treleani, Phys. Rev. D60 (1999) 054023. 
[10] FA. Ceccopieri, Phys. Lett. B697 (201 1) 482. 
[11] W. Schafer and A. Szczurek, Phys. Rev. D85 (2012) 094029. 
[12] A. Donnachie and PV. Landshoff, Phys. Lett. B296 (1992) 227. 
[13] M. Luszczak, R. Maciula and A. Szczurek, Phys. Rev. D85 (2012) 094034. 



6 



